Weak Measurement of Qubit Oscillations with Strong Response Detectors: Violation 
of the Fundamental Bound Imposed on Linear Detectors 
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We investigate the continuous weak measurement of a solid-state qubit by single electron transis- 
tors in nonlinear response regime. It is found that the signal-to-noise ratio can violate the universal 
upper bound imposed quantum mechanically on any linear response detectors. We understand the 
violation by means of the cross-correlation of the detector currents. 



Single electron transistor (SET) is a sensitive charge- 
state detector [H, H, ll], which promises the use for fast 
qubit read-out in solid-state quantum computation. For 
single-shot measurement, i.e., in one run the qubit state 
is unambiguously determined, an important figure of 
merit is the detector's efficiency, defined as the ratio of 
information gained time and the measurement induced 
dephasing time 0. In the weakly responding regime, it 
was found that SET has rather poor quantum efficiency 
[3, 0, Q. However, recent study showed that, for strong 
response SET, the quantum limit of an ideal detector can 
be reached, resulting in an almost pure conditioned state 



Rather than the single-shot measurement, a more im- 
plementable approach in experiment is the continuous 
weak measurement. This type of measurement allows the 
ensemble average of detector and qubit states, and the 
qubit coherent oscillation is read out from the spectral 
density of the detector. In continuous weak measure- 
ment, a remarkable result is the Korotkov-Averin (K- 
A) bound, originally stated as follows [3]: The interplay 
between the information acquisition and the backaction 
dephasing of the oscillations by the detector imposes a 
fundamental limit, equal to four, on the signal-to-noise 
ratio of the measurement. The limit is universal, e.g., 
independent of the coupling strength between the detector 
and system, and results from the tendency of quantum 
measurement to localize the system in one of the mea- 
sured eigenstates . In order to overcome the K-A boimd, 
special techniques such as the quantum nondemolition 
(QND) measurement Q, the quantum feedback control 
0, and the measurement with two detectors have 
been proposed. In this work we investigate the contin- 
uous weak measurement by strongly responding SETs 
d, [ll| . Remarkably, we find that for both models stud- 
ied in Refs. and [ill the signal-to-noise ratio (SNR) can 
violate the universal Korotkov-Averin bound. Interpre- 
tation and implication of this new result will be also pro- 
vided. 
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FIG. 1: Schematic model for a solid-state qubit measurement 
by SET. Model (I): the SET dot level is within the bias voltage 
for qubit state |6), but outside of it for state \a). Model (II): 
the SET dot level is between the Fermi levels for either \b} 
or \a), but with different coupling strengths to the leads, i.e., 
Fi/ii for \b), and F'^/^ for \a). 



Model. — Consider a charge qubit, say, an electron in 
a pair of coupled quantum dots, measured by a single 
electron transistor, as schematically shown in Fig. 1. The 
entire system is described by the following Hamiltonian 



H = Ho+H' 
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For simplicity, we assumed spinless electrons. The sys- 
tem Hamiltonian, Hs, contains qubit and the SET cen- 
tral dot, and their Coulomb interaction (the [/-term). 
For qubit, we assumed that each dot has only one bound 
state, i.e., the logic states |a) and \b) with energies Ea and 
Eh, and with a coupling amplitude fl. is the number 
operator of qubit state |a), which equals 1 for \a) occu- 
pied and otherwise. For the SET, al{ac) and d^^f^{dak) 
are the electron creation (annihilation) operators of the 
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central dot and reservoirs, ric = ajac is introduced as 
the number operator of the SET dot. Similar to previ- 
ous work, we assumed that the SET works in the strong 
Coulomb-blockade regime, with only a single level Ec in- 
volved in the measurement process. Finally, H' describes 
the tunnel coupling of the SET dot to the leads, with am- 
plitudes flxk- 

In this work, we consider two SET models as schemat- 
ically shown in Fig. 1. In model (I), which was studied 
in Ref. [Ill, the SET dot level is within the bias voltage 
if the qubit is in state \b), but it locates above the Fermi 
levels when the qubit state is switched to \a). For state 
1 6), a nonzero current lb flows through the SET; how- 
ever, for state |a), the SET current la is zero. Then, 
the qubit state can be discriminated from these different 
currents. In this model, the signal cmTcnt AI = \Ib — Ia\ 
is twice the average current / = {If, + la)/ 2. In this 
sense, it is not a weak response detector. In model (II), 
which allows to illustrate the crossover from weak to 
strong responses, the SET dot level is always between 
the Fermi levels of the two leads, for qubit either in state 
1 6) or in state |a), but with different coupling strengths 
to the leads, i.e., ^l{R} and ri(jj)- For the convenience 
of description, we further parameterize the tunnel cou- 
plings as TLir^) = (1 ± e)f i, r«(r'^) = (1 ± C)f fl, and 
7 = ffl/fi. Here, T^r^^ = {T^n^ + r^(^))/2 denote 
the average couplings, while ^ and C characterize the re- 
sponse strength of the detector to qubit. In this context, 
we would like to stress that most previous works were 
largely restricted in the weak response regime by assum- 
ing ^ ^ 1 and C ^ Ij except in Ref. |6j the quantum 
efficiency was investigated in strong response regime us- 
ing this model. 

Method. — In continuous weak measurement, the de- 
tector's output is characterized by the current and noise 
spectral density. For their calculation, the most efficient 
approach is the particle-number resolved master equation 
[2]. In obtaining it, the qubit and SET dot are regarded 
as the system of interest, while the two leads of the SET 
as an environment; and the tunnel coupling H' of SET 
is treated perturbatively as an interaction between them. 
Up to the dominant second-order of H' , following Ref.[T3 
we have 

- + Ai-Jp^---'^al] + H.c.}. 

(2) 

pina) ^Yie reduced density operator of the system con- 
ditioned on the electron number "nfl" tunnelled through 
the right junction (similar equation holds also for the 
left junction). For simplicity, throughout this paper 
we use the convention h = e ~ kg ^ I. In Eq. ^ 
the Liouvillian C is defined by £(• • • ) = [Hs, ■ • • ], and 
the operators ^^a^ = C^\±C)ac. The superoper- 
oters C^\±C) are the generalized spectral functions: 



C[^\±C) = J^^ dtc[^\t)e'^'^*, where the bath corre- 
lation functions C['^Ht) = {fl)^{t)fc\)B and &y^\t) = 

{fc\{t)flx)B, and the average {■■■)b = Ttb[{- ■ ■)pb], 
with pb the local thermal equilibrium state of the SET 
leads determined by the respective chemical potentials. 

Rich information is contained in the above particle- 
number resolved master equation, since the conditional 
density matrix p^'^'^'>{t) is directly related to the distribu- 
tion function, P{nii,t) = Tr[p^""' (i)], where the trace is 
over the system states. For instance, by virtue of this re- 
lation, the measurement current can be obtained as p^ : 

In{t) = Tr{n«p("«)(i)} = RcTr[(at^(^) - A[+^a\)p{t)l 

(3) 

where p{t) = X^hr /^'"^''(^) satisfies the usual uncondi- 
tional master equation, by summing the above Eq. (P| 
over nn. 

In continuous weak measurement, the detector's power 
spectral density contains very useful information of 
qubit's coherent oscillation. Formally, the noise spec- 
trum of current consists of three terms S{uj) = 
aSLito) + 0Sfi,{Lu) — a/3ci;^S'Ar(cj), with Sl/ii(uj) the noise 
of the left (right) junction current lL/R{t), and Sn{^^) 
the fluctuations of the electron number N{t) on the cen- 
tral dot of SET. a and /3 are two coefficients determined 
by the junction capacitances p^ . and satisfy a + (3 = 1. 
Further, for Sl/r{lo), it follows the MacDonald's formula 

Sxii^)^2cjJ^ dtsinujt-[{nl{t)) - {It)\ (4) 

where / is the stationary current and {n\{t)) = 
S„^n|Trp("^^ (<) = E„^n|P(n;)j, t). In practice, instead 
of directly solving P{nx,t), the reduced quantity {n\{t)) 
can be obtained more easily by constructing its equa- 
tion of motion [l^ , based on the particle-number resolved 
master equation Eq. Thus, we can obtain 

|K(i)) = Tr[2ji-)Q,(t) + 4+'^p% (5) 

where the particle-number matrix reads Qx{t) = 
n^^yO^"^^ (<) and p'** the stationary density matrix. 
Here the superoperator means 

4^\---)-l[4-x\---)al±ali...)A[X' 

+ a,(...)4A^t±4|)^(---)«c]. (6) 

Inserting Eq.(5) into Eq.(4) and transforming into the 
frequency domain, we obtain 

Sx{uj) = 4u;lm{TT[J^-'^qxiuj)]} + 2TrJ^+V'* - SnPS{uj), 

(7) 

in which Qx{ijj) = /o°° QA(i)e"^*. we can easily calcu- 
late Qxi^^) by solving a set of algebraic equations after 



3 



Fourier transforming the equation of motion of Q\{t), 
as have been clearly described in ReflTsl. For Sn{lu), 
which is the Fourier transform of the correlation function 
{{N (t) , N {0)}) , following Ref. [H the quantum regres- 
sion theorem gives SN{i^) = 2ReTr{N[(T(cj) + (t(— w)]}. 
ct(w) is introduced as the Laplace transform of a{t) = 
Tr B[U{t)Np''*PBW{t)], where U{t) = e-'"--^*. Obvi- 
ously, a{t) satisfies the same equation of the reduced 
density matrix p(t). The only difference is the initial 
condition, for a{t) which is a{0) = Np"*. 

Results. — For both models in Fig. 1, the states in- 
volved arc |1) = |0a), |2) = |06}, |3) = |la}, and 
|4) = In this notation |0(l)a(6)) means that the 

SET dot is empty (occupied) and the qubit is in state 
\a{b)). Applying Eq. ^ to model (I) results in 



p[r ^ =^n[p[T^ - piT^] + T.pir ^ + r^pir "'^ (8a) 

.(n 

P22 



P21 ^ P12 J ^ ^ LP22 + ^ RPii V^^) 
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^^n[piT^ - pi'-) 



^ RP3i 



{TR + TL)p'i^ 



(8c) 
(8d) 



r^p'r'-rKpir) (8e) 



P^^,"■^--^ie + u)piT^ + ^n[pir^-piTh 



L Ann) 
o P12 



O 'Psi 



(8f) 



Here, e = Ea-E^ and F^/fl = 2Tr\nL/ji\'^gL/^]i, with g^/^ 
the density of states of the SET leads. For simplicity, the 
assumption of wide-band limit implies il^L/R = ^L/B.ki 
and makes Tj^/jf energy independent. Also, low temper- 
ature and U Q were assumed to further simplify the 
equations. Similarly, for model (II), we have 

pjr' -mp[t^ - pir^] - r',p[r^ + r^nPiT''^ m 

(9b) 



P22" 

P12 ~ 'tpj2 



pi-Vr^P^r^ + F^pi: 



-P34 



2 2 

('ih)i , -p' „("h) p/ „("«) 



P43 J + ^ lPii ^ RP33 



p(r^=-z(6+c/)p^r'+*f^[pir'-pir^] 

,'£±±SAJnR) + „(nH) 
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(9c) 

(9d) 
(9e) 



(9f) 



Except for the conditions leading to model (II), other 
parameters are the same as above. 

In continuous weak measurement of qubit oscillation, 
the signal is manifested as a peak in the noise spectrum 
at the qubit oscillation frequency 2fl, while the measure- 
ment effectiveness is characterized by the signal-to-noise 
ratio (SNR). i.e., the peak-to-pedestal ratio. We denote 




FIG. 2: Signal-to-noise ratio: (A) for model (I), and (B) for 
model (II). For model (I), we take Fl = F as the energy unit, 
and assume that Plir) = ±50r, n = 2r and U = SOP. For 
model (II), we use Fl = F as the energy unit, and assume 
that Q = f, U = 50r, Tr = 30f , and a = /3 = 1/2. Also, 
zero temperature and Ea = Et are assumed. 



the noise pedestal by S'p, and obtain it conventionally 
from S{uj 00). In Fig. 2 we show the dependence of 
the SNR on detector's configuration symmetries. 

The result of model (I) is shown in Fig. 2(A), where we 
see that both the tunnel- and capacitive-coupling symme- 
tries crucially affect the measurement effectiveness. For 
the effect of tunnel coupling asymmetry Th/Tl, the basic 
reason is that, with the increase of Fj^/Fi, the interaction 
time of the detector electron with the qubit is decreased. 
Thus the detector's back-action is reduced and the SNR 
is enhanced For the effect of capacitive coupling, 

its degree of asymmetry affects the contribution weight 
of the cross-correlation between lL{t) and lR{t) to the 
entire circuit noise. Specifically, the cross-correlation has 
more important contribution for more symmetric cou- 
pling, as shown in Fig. 2(A) by the a-dependencc. This 
is because, as we shall demonstrate below, the cross- 
correlation has much higher peak-to-pedestal ratio than 
the auto-correlation. 

An unexpected feature observed in Fig. 2(A) is that 
under proper conditions, say, the symmetric capacitive 
coupling and strongly asymmetric tunnel coupling, the 
SNR can exceed "4" , which is the upper bound quantum 
mechanically limited on any linear response detectors 0]. 
However, to our knowledge, whether this upper bound is 
applicable to nonlinear response detector is so far unclear 
in priori, since in this case the linear response relation 
between current and qubit state breaks down, then the 
subsequent Cauchy-Schwartz-incquality based argument 
leading to the upper bound "4" docs not work [10l |. 

To support the above reasoning, we further check 
model (II). The result is presented in Fig. 2(B). As ex- 
plained in the model description, the parameters ^ and 
C used here characterize, respectively, the left and right 
tunnel-coupling responses to the qubit states. Shown in 
Fig. 2(B) is for an asymmetric tunnel coupling detector, 
with 7 = Th/Tl = 30, which can lead to higher SNR, 
because of the weaker back-action from the detector, sim- 
ilar to model (I). Here we find that the SNR is insensitive 
to the right junction response C, but sensitive to the left 
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one ^. Again, in this model, we observe that the SNR 
can violate the K-A bound "4" in the strong response 
regime. 




6 12 18 24 300_ 6 Jl2 18 24 30 
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FIG. 3: Signal-to-noise ratio versus tunnel-coupling asymme- 
try, Vr/Vl for model (I), and Vr/Vl for model (II). In (A) 
the solid and dashed lines are the result in the presence of 
cross correlation, while the dotted and dot-dashed lines are 
the result after removing it. In (B) the mere cross correlation 
is plotted. Sp is the pedestal noise of the entire circuit cur- 
rent. ^ = C = 0.9, other parameters are the same as in Fig. 
2. 



Understanding the Violation of the K-A Bound. — 
Since I{t) — alhit) + piR{t), the current corre- 
lator (/(t)/(0)) contains the component SLR{t) = 
{lL{t)lR{0) + //f(t)/L(0)}, i.e., the cross-correlation. 
Also, in the previous results, we see that for more sym- 
metric capacitive coupling the SNR is larger, and reaches 
the maximum at a = (3 = 1/2. This feature indicates 
that the cross-correlation has an effect of enhancing the 
SNR. Indeed, for the SET detector, both the left and 
right junction currents {II and Ir) contain the informa- 
tion of qubit state, so their "signal" parts are correlated. 
This leads to a heuristic opinion that views the two junc- 
tions as two detectors, like the scheme of qubit measure- 
ment by two point-contacts proposed recently by Jordan 
and Biittiker , where they found that the SNR of the 
cross correlation can strongly violate the K-A bound, be- 
cause of the negligibly small pedestal of the cross noise. 
In our case, since Ihit) and lR{t) are subject to a con- 
straint from charge conservation, the cross noise back- 
ground of Ihit) and lR{t) does not vanish in p rinciple, 
unlike the two independent QPC detectors Never- 
theless, the pedestal of the cross noise of SET is much 
smaller than that of the auto-correlation, which leads to 
an enhanced SNR in the spectral density of the total cir- 
cuit current, and to the violation of the K-A bound, as 
clearly shown in Fig. 3(A). For comparative purpose, in 
Fig. 3(B) we plot the SNR of the cross correlation, scaled 



by the noise pedestal Sp of the circuit current. 

In Fig. 4 the spectral density of the cross correlation, 
scaled by its own noise pedestal, is shown representa- 
tively. As mentioned above, since at high frequency limit 
the cross noise pedestal is negligibly small, here we ar- 
tificially (but more physically in some sense) define the 
pedestal at a finite frequency, e.g., twice the qubit oscil- 
lation frequency. Obviously, the giant SNR of the cross- 
correlation has drastically violated the K-A bound. This 
result indicates that in qubit measurement by SET one 
can exploit the cross correlation, rather than the auto one 
as usual, to probe the coherent oscillations. In practice, 
such scheme is simpler than the technique of QND mea- 
surement Q , and holds the most advantages of SET over 
QPC. In recent years, the cross-correlation in mesoscopic 
transport is an extensive research subject. Its measure- 
ment in experiment is also possible, for instance, using 
the nearby-QPC counting technique p^ . 
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FIG. 4: Spectral density of the cross-correlation scaled by 
its own pedestal, here which is defined at twice the Rabi fre- 
quency of the qubit oscillations. Parameters for model (I) 
in (A): Fl = O.OSr, Fr = O.SF, and = 5r. Parameters 
for model (II) in (B): Fl = O.OSF, Tr = O.SF, Q = 5F, and 
^ = ( — 0.9. Fin this figure is used as an energy unit, other 
conventions are the same as in Fig. 2. 



Summary. — We have investigated the continuous weak 
measurement of qubit oscillations by nonlinear response 
SET, and demonstrated that the signal-to-noise ratio can 
violate the universal Korotkov-Avcrin bound. The vi- 
olation has been understood by the role of the cross- 
correlation of the detector's currents. This interesting 
interpretation also leads to useful implication to experi- 
ment. 

Acknowledgments. This work was supported by the 
National Natural Science Foundation of China under 
grants No. 60425412 and No. 90503013, the Major State 
Basic Research Project under grant No.2006CB921201. 



5 



[1] M. H. Devoret and R. J. Schoelkopf, Nature (London) 

406, 1039 (2000). 
[2] Yu. Makhlin, G. Schon, and A. Shnirman, Rev. Mod. 

Phys. 73, 357 (2001). 
[3] W. Lu, Z. Ji, L. Pfeifer, K.W. West, and A.J. Rimberg, 

Nature 423, 422 (2003). 
[4] A. N. Korotkov, Phys. Rev. B 63, 085312 (2001); ibid 

63, 115403(2001). 
[5] D. Mozyrsky, I. Martin, and M. B. Hastings, Phys. Rev. 

Lett. 92, 018303 (2004). 
[6] N. P. Oxtoby, H. M. Wiseman, and H. B. Sun, Phys. Rev. 

B 74, 045328 (2006). 
[7] A. N. Korotkov and D. V. Averin, Phys. Rev. B 64, 

165310 (2001). 

[8] D. V. Averin, Phys. Rev. Lett. 88, 207901 (2002); A. N. 
Jordan and M. Biittiker, Phys. Rev. B 71, 125333 (2005). 



[9] S. K. Wang, J. S. Jin, and X. Q. Li, Phys. Rev. B 75, 
155304(2007) 

[10] A. N. Jordan and M. Biittiker, Phys. Rev. Lett. 95, 
220401 (2005). 

[11] S. A. Gurvitz and G.P. Berman, Phys. Rev. B 72, 073303 
(2005). 

[12] X. Q. Li, P. Cui, and Y. J. Yan, Phys. Rev. Lett. 94, 
066803 (2005). 

[13] J. Y. Luo, X. Q. Li, and Y. J. Yan, Phys. Rev. B 76, 
085325 (2007). 

[14] S. Gustavsson, R. Leturcq, B. Simovi, R. Schleser, T. 
Ihn, P. Studerus, K. Ensshn, D.C. Driscoll, and A.C. 
Gossard, Phys. Rev. Lett. 96, 076605 (2006); T. Fuji- 
sawa, T. Hayashi, R. Tomita, and Y. Hirayama, Science 
312, 1634 (2006). 



